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ABSTRACT 

ERS-1 SAR imagery of lake ice growing on shallow tundra lakes in 
northern Alaska shows interesting radar backscatter variations. Based on the 
analysis of ice cores from these lakes, a multi-layer backscatter model 
comprised of tlie following elements has been developed: I) specular air-ice. 
ice-water and ice-frozen soil boundaries; 2) an ice layer of variable thickness; 
3) ice sub-layers with air inclusions of variable density, size and shape, 
including spheres, prolate spheroids, and cylinders of finite length. 
Preliminary model results confirm that backscatter is a sensitive function of 
the presence of a specular ice-water interface, with a roughly 40-times 
greater reflectivity than from an ice-frozen soil interface. The model has also 
been tested using bubble data derived from ice cores in April 1992. The 
modelled backscatter is compared with backscatter derived from ERS-I SAR 
images obtained at the same time as the fieldwork. 
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I. INTRODUCTION 

The North Slope of Alaska is a large expanse of low-lying tundra with 
many shallow lakes covering more than 40% of the surface area. During the 
1970s, in late winter, SLAR images of these lakes showed interesting 
variations of backscatter intensity, with areas of low backscatter at lake 
margins believed to represent ice that was frozen to the lake bed, while areas 
of high backscatter were considered to represent floating ice that contained 
tubular bubbles which acted as forward scatterers |Sellman et al.,1975; 
Elachi et al.,1976: Weeks et al , 1 977. 1978. 198 1 1. 

From September 1991 to April 1992, backscatter intensity variations from 
shallow tundra lakes near Barrow, NW Alaska, were studied using C-band 
SAR data from the ESA Remote-Sensing Satellite-1 (ERS-I). The SAR 
data were processed at the Alaska SAR Facility (ASF) and backscattering 
coefficients were derived for a number of lakes. Backscatter intensity varied 
from -19 dB to -ft dB. Held measurements in April 1992 confirmed that the 
highest values were associated with floating ice containing tubular bubbles, 
while the lowest backscatter values were associated with ice frozen to the 
lake bed. but which also contained tubular bubbles (Jeffries et al., 1 993 (. 

This paper repons on the development of the first numerical model of C- 
band backscatter from such ice covered, shallow tundra lakes. Backscatter is 
simulated for grounded and floating ice of variable thickness and 
stratigraphic characteristics. Also, using ice physical characteristics observed 
and measured in late April as inputs to the model, simulated backscatter is 
compared with backscatter derived from ERS-I SAR images obtained al a 
similar time. 


2. MODELING 


2.1 Multi-layer consideration 

In April 1992, field measurements were made at 1 1 locations near Barrow, 
NW Alaska. Ice cores were also obtained from 10 sites and taken to the 
laboratory for detailed ice core analysis, which indicated the presence of three 
ice types: 1) granular ice containing roughly spherical bubbles, with radii 
much smaller than a wavelength (« 5.7 cm); 2) clear bubble-free ice; and, 
3) ice containing tubular bubbles resembling thin cylinders. The latter, 
oriented vertical to the ice surface, also had radii much smaller than a 
wavelength, but with lengths ranging from 1.5 cm and 9.1 cm. The number 
of tubular bubbles per square meter, i.e. bubble density or ice porosity, was 
quite variable (Jeffries et al., 1993|. 

The typical stratigraphy of the lake ice is shown in Ftg. 1. The clear ice 
layer and one or more tubular bubble layers were found at all sites, while the 
granular ice layers- were found at 6 of the 10 sites. Each layer had a different 
thickness, bubble density and bubble size, especially in the tubular bubble 


layers where the bubble density and bubble size varied considerably with 

depth. The snowcover was of variable depth (0-80 mm) and dry. 

Considering the features described above, our basic lake ice model is 

comprised of the following elements: 

(1) three main layers (either air-ice-water or air-ice- frozen soil); 

(2) sub-layers in the ice layer with each sub-layer having a different bubble 
density and bubble size (a bubble density of zero corresponds to clear 
ice); 

(3) the boundaries between the main layers are considered to be specular; 

(4) reflections occur at the main layer boundaries and at the air bubbles; and. 

(5) dielectric parameters of the ice layers change with the bulk ice density, 
i.e.. with the density of bubbles present in the ice. 


Because the surface roughness on these lakes is generally considered to be 
small compared to other sources of scattering, surface scattering is neglected. 
The dry snowcover is also neglected. For dry snow, the dielectric mismatch 
at the air-snow boundary is very small and the reflection at the air-snow 
boundary can be ignored. In addition, since the snowcover was thin and 
snow particles very small, the contribution of the snowcover to the 
backscattering coefficient may be neglected at both L- and C-band 
frequencies |U!aby et al.. I986|. It should be noted that reflections are added 
incoherently in the model, because the length of each path is widely 
distributed. Calculations are made on the basis of power. 
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The scattering coefficient of an extended target in a given direction is the 
ratio of the total scattered power from an equivalent isotropic scatterer which 
generates the same scattered power density in the direction to the total power 
of the illumination area. If the incident and scattering intensities are given, 
the backscattering coefficient would be as follows: 

f 

where 0 O is the incident and scattered angle, /' is the scattered p-polarized 

intensity, and l‘ q is the incident q-polarized intensity. When reflections by 

bubbles in the ice occur, the scattered power is propagated to the 
backscattered direction, as shown schematically in 
Fig. 2. In the general case, the bubble density and bubble size vary 
according to the depth of ice. Thus, it is assumed that the ice is composed of 
n-layers and each layer has a different bubble density and bubble size. The 
downward and the upward intensities (£>,0 ) at boundaries of layer i are 
defined as shown in Fig. 3. The relation of these intensities is given by 


( 2 ) 


where 

D. = [D* D-], D,_, =[0*. £>,:,]*. u, =[£/; u;],i 

Here t indicates a transposed matrix, and M d , M b , M c and M d are 2 x 2 
matrices, of which the elements are given by the calculation for the first 
order scattering. For convenience in the multi-layer treatment, the following 
matrix for each layer is calculated. 
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In order to obtain the relation between the incident and scattering intensities 
at the air-ice boundary, this matrix for each layer is derived and multiplied, 
and finally, the matrix for entire ice layer K derived according to. 
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where '* ,lle element of the matrix in row I and column m. Because 
there is no upward intensity from the bottom layer (water or frozen soil ) 
and no left downward intensity in the upper layer (air), the following 
conditions must be satisfied, 

d; = i‘.d;=u; =u;=o 


Using the above condition, the ratio of backscattering intensity to the 
incident intensity is given by. 
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Substituting Eq.(5) to Eq.(l), we obtain the backscattering coefficient as 
follows. 
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2.2 Bubble scattering 

In order to obtain each element in the 4 X 4 matrix of Eq.(2), the reflectivity 
of each boundary, the volume extinction coefficient, and the bistatic cross 
section for bubbles are required. The power reflectivity can be calculated by 
applying Fresnel's law using the average permittivity of the ice layer, which 
is considered to be a “two phase mixture" composed of pure ice and air 
bubbles IdeLoor. 1 968 1. The volume extinction coefficient is the sum of the 
volume absorption coefficient and volume scattering coefficient. 


grounded ice is -10 to -12 dB [Jeffries et al„ 1 993 1. Since the power 
reflectivity for an ice-water interface is approximately 40 limes that of an 
ice-soil interface, the simulation confirms that backscatter is significantly 
changed by the presence or absence of a specular ice-water interface. 

Backscatter change as a function of ice thickness for an incidence angle of 
23 is shown in Fig.6. In the simulation, the ice had 10cm thick granular 
layer, a 60cm thick clear layer and 20-200cm thick tubular bubble layer. The 
porosity and bubble size were the same as those used in the previous 
simulation. The increase in backscatter with increasing ice thickness, 
ranging from -14 dB to -6 dB, is similar to that observed by SAR during the 
period from October 1991 to April 1992 [Jeffries et al„ 1993). 

Other simulations showed that backscatter increases as bubble radius, bubble 
length and bubble density increase. 

4. COMPARISON BETWEEN ERS-1 SAR AND 
MODELLED BACKSCATTER 

The final part of this paper compares simulated and observed backscatter. 
For the simulation, the ice was divided into sub-layers, and the porosity, 
mean bubble radius and mean bubble length in each sub-layer as measured in 
ice cores obtained in April 1992 were used as model input. In order to obtain 
the observed data, the digital number (DN) for a specific area of a lake was 
selected from low-resolution data processed at ASF and converted into 
backscattering coefficients. 

The result of the comparison between SAR-derived backscatter and the 
model output is shown in Fig. 7. At most sites, the model underestimates 
the backscatter, particularly at sites 7, 8 and 1 1. The correlation coefficient 
for a linear regression (Fig. 7) of the two data sets is 0.73 and the model has 
-5.5 dB offset relative to the SAR data. 


The air bubbles in the granular ice layer are modelled as spheres with a size 
much smaller than a wavelength. In the tubular bubble ice layer, the air 
bubbles are modelled as either prolate spheroids or cylinders of finite length 
according to the length of bubbles. Because the radii of the bubbles are 
much smaller than a wavelength for all cases, low frequency approximation 
is applied to obtain the bistatic cross section. The scattering formulation 
used in this model is that given by Ruck et al. (1970). 


The geometry of the problem is shown in Fig. 4. The bistatic VV cross 
sections for spheres cf'Vv/,, y/,,0), for prolate spheroids <f''(i//,, 
and for cylinders of finite length o*(. g>„ ys,,<p) are given as follows; 
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where k is the wave number, a is the bubble radius, E is the relative 
permittivity, h is the half length of bubble, a, and b t are the semi-axes of 
the equivalent spheroid, /„ and / 6 are elliptic functions, and <f( I/. , 0) is 
the scattering width of an infinitely long cylinder at oblique incidence. 

3. SIMULATION 

Backscatter from both floating ice and grounded ice have been calculated 
using this model. For the simulation, the following ice thickness and 
stratigraphic characteristics were used. The grounded ice consisted of a 10cm 
thick granular layer, a 15cm thick clear layer and a 30cm thick tubular 
bubble layer. The floating ice consisted of a 10cm thick granular layer, a 
60cm thick clear layer and a 1 40cm thick tubular bubble layer. In the 
granular layer the porosity was I % and the bubble diameter was 2mm, while 
in the tubular layer the porosity was 3% and the bubbles were 1mm in 
diameter and 5cm long. In the simulations the real and imaginary C-band 
relative dielectric constants for the different materials were as follows; Air 
1.0. 0; Ice 3.15, 9.00E-4; Water 65.00, 35.00; Soil 5.00, 0.50. 

The results of the simulation are plotted in Fig. 5 which clearly shows the 
large difference in backscatter intensity between floating and grounded ice. 
At the ERS-1 incidence angle (23°) the backscatter difference is about 7.5 
dB. Using ERS-1 SAR data the backscatter difference between floating and 


S. SUMMARY 

On the basis of ice characteristics obtained from ice core measurements, a 
numerical model for C-band backscatter from ice growing on shallow tundra 
lakes has been developed. The model explains the backscatter difference 
between grounded ice and floating ice, and also the observed increase in 
backscatter as the ice grows thicker during the winter. Backscattering 
coefficients derived from ERS-1 SAR data were greater than those estimated 
by the model. The reasons why the model under-estimates backscatter 
remain to be determined. Further model improvement, such as incorporating 
a surface roughness model and multi-order bubble scattering, as well as 
additional field measurement are needed. 
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Fig. I. Typical structure of lake ice as observed in Northern Alaska 
f P P 



80 :incident and scattered angle 
d : ice thickness 
I' .-incident intensity 
1* : scattered intensity 

Ftg. 2. Schematic figure of lake ice reflection including direct backscattering 
and forward scattering from bubbles. 



Fig. 5. Backscattering coefficient at C-band (VV) versus angle of incidence 
for the floating ice and the grounded ice. 
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Ftg. 3. Multi-layer ice structure 
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Fig. 6. Backscattering coefficient at C-band (VV) versus the thickness of an 
ice layer with tubular bubbles. 



Fig. 7. The backscattering coefficient observed by ERS-I SAR versus the 
backscattering coefficient estimated by lake ice model. Numbers 
displayed on the figure indicate the site number. The dotted line is the 
result of linear curve fitting. 


Fig. 4. Bubble geometry 
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